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The influence of short-range correlations on the momentum and energy distribution of nucleons 

^ , in nuclei is evaluated assuming a realistic meson-exchange potential for the nucleon-nucleon inter- 

O ' action. Using the Green-function approach the calculations are performed directly for the flnite 

' ^ I nucleus ^''O avoiding the local density approximation and its reference to studies of inflnite nuclear 

matter. The nucleon-nucleon correlations induced by the short-range and tensor components of the 

interaction yield an enhancement of the momentum distribution at high momenta as compared to 

the Hartree-Fock description. These high-momentum components should be observed mainly in 

►^ . nucleon knockout reactions like (e, e'p) leaving the flnal nucleus in a state of high excitation energy. 

\jr\ ' Our analysis also demonstrates that non-negligible contributions to the momentum distribution 

f^ , should be found in partial waves which are unoccupied in the simple shell-model. The treatment of 

t*^ ' correlations beyond the Brueckner-Hartree-Fock approximation also yields an improvement for the 

calculated ground-state properties. 
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I. INTRODUCTION 

Q 

^ \ Many properties of nuclei can be understood within the independent particle model (IPM). In the IPM the nucleus is 

considered to be a system of nucleons moving without residual interaction in a mean field or single-particle potential. 
The single-particle potential is either adjusted in a phenomenological way (assuming e.g. a Woods-Saxon shape) 
k^ , or evaluated from empirical effective interactions like the Skyrme forces within the Hartree-Fock approximation. 
^ ■ Attempts to employ realistic nucleon-nucleon (NN) interactions, which reproduce the NN scattering data, directly in 
^ \ such a scheme fail badly: typically one does not even obtain any binding energy in this approach. This result is due 
to the strong short-range and tensor components, which are typical for realistic interactions and induce corresponding 
NN correlations in the nuclear wave function, which cannot be described by the IPM or the Hartree-Fock approach. 
Various tools have been developed to account for these strong short-range correlations. These include variational 
calculations assuming Jastrow correlation functions M, the correlated basis function method (CBF) B], the "expo- 
nential S" method |^, the Brueckner-Hartree-Fock (BHF) approximation Q and the self-consistent Green function 
approach [||. 

Considerable effort has also been made to find a nuclear property which is experimentally accessible and reflects 
the effects of NN correlations in a clear manner. A candidate for such an observable is the momentum distribution of 
nucleons in a nucleus. This momentum distribution can be written as 

= E/.,,.(2j + 1) {HWu^rakl.r \H) ■ (1) 

Here h&o ) represents the ground state of the nucleus under consideration (with A nucleons) and aj.^ ■ (akijr) denotes 
the creation (annihilation) operator for a nucleon with orbital angular momentum /, total angular momentum j, 
isospin r and momentum k. The momentum distributions for the partial waves, nijr{k), in Eq. (n^) can be rewritten 
by inserting a complete set of eigenstates |\E'^^^) for the system with A— \ nucleons 
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nyrW=^|«-l|aM,r|*^)f. (2) 
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In the IPM the sum in this equation is typically reduced to one term, if (Z, j, r) refer to a single-particle orbit occupied 
in |^(?)- Eq. (||) then yields the square of the momentum-space wave function for this single-particle state. The 
contribution nijr{k) vanishes in the IPM if no state with quantum numbers {ItJ,t) is occupied. If correlations are 
present beyond the IPM approach this simple picture is no longer true and the determination of the momentum 
distribution n{k) requires both in experimental as well as theoretical studies the complete knowledge of the nucleon- 
hole spectral function 

S^,Ak,E) = J2\{^t'\aH,r \^^)\'S{E - {E^ ~ Et')) (3) 

for all energies E and all sets of discrete quantum numbers (Z, j, r). Note that the energy variable E in this definition 
of the spectral function refers to the negative excitation energy of state n in the A — 1 system with respect to the 
ground-state energy of the nucleus with A nucleons {Eq). The spectral function is experimentally accessible by 
analyzing nuclcon knockout experiments like (e,e'p). The momentum distribution nijT-{k) is obtained by integrating 
the spectral function over energies E from — c» to the Fermi energy ep = Eq — Eq^ with Eq ~ denoting the energy 
of the ground state for A— 1 nucleons. One important aim of our studies is to investigate how short-range correlations 
modify the spectral function at various energies as compared to the IPM. For example, can one expect to observe 
effects of short-range correlations in knockout experiments with small energy transfer? 

Microscopic calculations of the spectral function and the momentum distribution based on realistic nuclear hamil- 
tonians have mainly been performed for very light nuclei (A < 4) &-Pl or nuclear matter pO[-p7|. Results for heavier 
nuclei are typically derived from investigations of nuclear matter assuming a local density approximation [[l8|-pi[. 
Recent variational calculations for ^^O yield the momentum distribution 123] and the p-shell quasihole wave functions 
p3[ but not the complete energy dependence of the hole spectral functions. 

The present investigation determines the spectral function and the corresponding momentum distribution directly 
for finite nuclei without employing the local density approximation. Calculations are performed for the nucleus ^^O 
assuming a realistic meson-exchange potential p4] for the NN interaction. The spectral function is derived from the 
Lehmann representation of the single-particle Green function. This Green function solves the Dyson equation with a 
self-energy calculated by techniques as described in ref. [|5|. A few results concerning the spectral function for the 
Pi/2 partial wave have been discussed already in a brief report [2^ . 

Special attention is paid to the effect of correlations on the spectral function at different energies. We find that 
clear indications of the short-range NN correlations are obtained by studying the spectral function at very negative 
energies, which in nucleon knockout experiments correspond to excitation energies of around 100 MeV and more in 
the remaining nucleus. The resulting Green function is also used to study the effects of correlations beyond the BHF 
approach on the binding energy and radius of the nuclear ground state. 

After this short introduction we describe the techniques used to evaluate the spectral functions and momentum 
distributions in section 2. The results of our numerical studies are presented in section 3 and the final section 4 
summarizes the main conclusions of this work. 

II. EVALUATION OF THE SPECTRAL FUNCTIONS 



The spectral function for the various partial waves, Sijr{k^E), (see Eq. (0)) can be obtained from the imaginary 
part of the corresponding single-particle Green function or propagator gij{k, k; E). Note that here and in the following 
we have dropped the isospin quantum number r. Ignoring the Coulomb interaction between the protons the Green 
functions are identical for N — Z nuclei and therefore independent of the quantum number r. The single-particle 
propagator can be obtained by solving the Dyson equation 

gij{ki,k2;E)^gl^\ki,k2;E)+ I dk3 j dk^g^f {ki,k:,-E)/\i:ij{k^M;E)gi,{ki,k2;E), (4) 

where g^^' refers to a Hartree-Fock propagator and ASjj represents contributions to the real and imaginary part of 
the irreducible self-energy, which go beyond the Hartree-Fock approximation of the nucleon self-energy used to derive 
g'^^\ The definition and evaluation of the Hartree-Fock contribution as well as the calculation of AE are presented in 
the next subsection. The methods used to solve the Dyson equation (0) and to extract spectral functions as well as 
momentum distributions are described in subsection 2.2. 



A. Nucleon Self-energy E 

The calculation of the self-energy is performed in terms of a G-matrix which is obtained as a solution of the 
Bethe-Goldstone equation for nuclear matter 

{k'l'SJsKLT\ G \k"l"SJsKLT) = {k'l'SJsKLT\ Vnn \k"l"SJsKLT) 

+ J2 [ ^^'^^ {k'l'SJsKLT\ Vnn \klSJsKLT) 
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{klSJsKLT\ G \k"l"SJsKLT) ^%P-^. (5) 

In this equation fc, fc', and k" denote the relative momentum, Z, V , and I" the orbital angular momentum for the 
relative motion, K and L are the corresponding quantum numbers for the center of mass motion, S and T denote 
the total spin and isospin of the interacting pair of nuclcons and by definition the angular momentum J5 is obtained 
from coupling the orbital angular momentum of relative motion and the s pin S. For the bare NN interaction Vnn we 
have chosen the One-Boson-Exchange potential B defined by Machleidt ( (24|, Tab. A. 2), m represents the mass of the 
nucleon, and the Pauli operator Q is approximated by the so-called angle-averaged approximation for nuclear matter 
with a Fermi momentum kp = lAfm~^. This roughly corresponds to the saturation density of nuclear matter. The 
starting energy ujnm has been chosen to be -10 MeV. The choices for the density of nuclear matter and the starting 
energy are rather arbitrary. It turns out, however, that the calculation of the Hartree-Fock term is not very sensitive 
to this choice [ p7| . Furthermore, we will correct this nuclear matter approximation by calculating the 2-particle 
1-hole (2plh) term displayed in Fig. |lp directly for the finite system, correcting the double-counting contained in the 
Hartree-Fock term (see discussion below) . 

Using vector bracket transformation coefficients ||28| , the G-matrix elements obtained from (g) can be transformed 
from the representation in coordinates of relative and center of mass momenta to the coordinates of single-particle 
momenta in the laboratory frame in which the 2-particle state would be described by quantum numbers such as 

\k1l1j1k2l2j2JT) , (6) 

where ki, li and ji refer to momentum and angular momenta of particle i whereas J and T define the total angular 
momentum and isospin of the two-particle state. It should be noted that Eq. (1q) represents an antisymmetrized 
2-particle state. Performing an integration over one of the fc^, one obtains a 2-particle state in a mixed representation 
of one particle in a bound harmonic oscillator while the other is in a plane wave state 
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\n1l1j1k2l2j2JT) = / dkikfRn^j^ 
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(aki) \k1l1j1k2l2j2JT) .. (7) 



Here Rm.ii stands for the radial oscillator function and the oscillator length a — 1.72 fm~^ has been selected. 
This choice for the oscillator length corresponds to an oscillator energy of hujosc — 14 MeV. Therefore the oscillator 
functions are quite appropriate to describe the wave functions of the bound single-particle states in ^®0. Using the 
nomenclature defined in Eqs. (ph - ffl) our Hartree-Fock approximation for the self-energy is easily obtained in the 
momentum representation 



Sgj;(fei,fci)= o^o.. ^1^ L {■2J+lK'2T+l){k,hjin2l2j2JT\G\k[hnn2l2j2JT). (8) 

n2l2J2JT 



The summation over the oscillator quantum numbers is restricted to the states occupied in the IPM of ^^O. This 
Hartree-Fock part of the self-energy is real and does not depend on the energy. 

The terms of lowest order in G which give rise to an imaginary part in the self-energy are represented by the 
diagrams displayed in Figs, pp) and^c), refering to intermediate 2plh and 2-hole 1-particle (2hlp) states respectively. 
The 2plh contribution to the imaginary part is given by 

Wlf^\k,,k[;E) = ^^^ Y.n2i2n T.IL EjjsST I k'dk J K^dK{2J + 1){2T + 1) 

X {kihjin2l2J2JT\ G \klSJsKLT) {klSJsKLT\ G \k[hj in2l2j2 JT) 

xn5[E + e,,,j,,,-£-^), (9) 

where the "experimental" single-particle energies f-n^hh ^^^ vsed for the hole states, while the energies of the particle 
states are given in terms of the kinetic energy only. The expression in Eq. (J^) still ignores the requirement that 



the intermediate particle states must be orthogonal to the hole states, which are occupied for the nucleus under 
consideration. The techniques to incorporate the orthogonalization of the intermediate plane wave states to the 
occupied hole states are discussed in detail by Borromeo et al. [g5[. The 2hlp contribution to the imaginary part 
Wi^-^{ki, k[;E) can be calculated in a similar way (see also [p5[). 

Our choice to assume pure kinetic energies for the particle states in calculating the imaginary parts of W'^^^^ (Eq. 
(y)) and W"^^^^ may not be very realistic for the excitation modes at low energy. Indeed a sizeable imaginary part 
in W'^'^'^P is obtained only for energies E below -40 MeV. As we are mainly interested, however, in the effects of 
short-range correlations, which lead to excitations of particle states with high momentum, the choice seems to be 
appropriate. A different approach would be required to treat the coupling to the very low-lying two-particle-one-hole 
and two-holc-onc-particle states in an adequate way. Attempts at such a treatment can be found in Refs. |29|-|32[|. 

The 2plh contribution to the real part of the self-energy can be calculated from the imaginary part VF^'' using a 
dispersion relation |3^ 

y^T(kuK;E) = f y_ '"V.!'^" ' dE'^ (10) 

where P means a principal value integral. From the (5-function in Eq. (H) one can see that W'^'^^^ is different from zero 
only for positive values of E' . Since the diagonal matrix elements are negative, the dispersion relation (|l^) implies 
that the diagonal elements of y^pi/i .^jjj i^g attractive for negative energies E. They will decrease and change sign 
only for large positive values for the energy of the interacting nucleon. 

Since the Hartree-Fock contribution E^^ has been calculated in terms of a nuclear matter G-matrix, it already 
contains 2plh terms of the kind displayed in Fig. ^). Therefore one would run into problems of doublecounting if 
one simply adds the real part \/2pih ^^ ^.j^g Hartree-Fock self-energy. Notice that Yi^^ does not contain any imaginary 
part because it is calculated with a nuclear matter G-matrix at a starting energy for which G is real. In order to 
avoid such an overcounting of the particle-particle ladder terms, we subtract from the real part of the self-energy a 
correction term, which just contains this contribution calculated in nuclear matter. This correction Vc is given by 

^^(^^' ^i) = 4(2j7+T) ^"-'^^^ ^'^ ^''^'^ ^ ^^'^^ ^ K'dK{2J + 1)(2T + 1) 

{kihjin2l2J2JT\ G \klSJKLT) {klSJKLT\ G \k[hjin2l2J2JT) 

^ ~ 5?^ — E^' UJ-J 
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with the same starting energy ujnm and the Pauli operator Q as used in the Bethe-Goldstonc equation (0). 
A dispersion relation similar to Eq. ( |lO| ) holds for ^2/iip ^^^ p^2hip 

P r°° W^''^P{ki,k[;E') 
V^^l^{,kuk[,E) = t: / _Jin±lL^^^dE', (12) 



^ 7-00 E'-E 

Since M^^^ip jg positive (at least its diagonal matrix elements) and different from zero for negative energies E' only, 
it is evident from Eq. (|l^) that ]/2pi/i jg repulsive for positive energies and decreases with increasing energy. Only for 
large negative energies it becomes attractive. 

Summing up the various contributions we obtain for the self-energy the following expressions 

E = E^^ + AE 

= E^^ -f (y2plft _Y^j^ y2hlp-^ _^ ^ {W'^plh _,_ ^r2hlp■^ (^3) 



B. Solution of the Dyson equation 

After we have determined the various contributions to the nucleon self-energy, we now want to solve the Dyson 
equation (H) for the single-particle propagator. In order to discretize the integrals in this equation we consider a 
complete basis within a spherical box of a radius i?box- This box radius should be larger than the radius of the 
nucleus considered. The calculated observables are independent of the choice of Rhox, if it is chosen to be around 15 
fm or larger. A complete and orthonormal set of regular basis functions within this box is given by 

$»ijm(r) = {r\kdjm) = N,iji{hr)yij„,{9<l)) (14) 



In this equation yijm represent the spherical harmonics including the spin degrees of freedom and ji denote the 
spherical Bessel functions for the discrete momenta hi which fulfill 



Using the normalization constants 



Ji(fc<^box) = 0. (15) 



■^ -, for/>0 



'RLJl-l{k^B■bo^) 



Nu = v: -— — ^ , ^ (16) 



T-yi 



for ; = 



the basis functions defined in Eq. (Ol) are orthogonal and normalized within the box 

flbox 

(fr {ki>l'j'm\v) {r\kjjm) = 5ii'5w5jj>5mm' ■ (17) 



Note that the basis functions defined for discrete values of the momentum fc^ within the box differ from the plane 
wave states defined in the continuum with the corresponding momentum just by the normalization constant, which is 
J^j-K for the latter. This enables us to determine the matrix elements of the nuclcon self-energy in the basis of Eq. 
(14) from the results presented in the preceeding section. 

As a first step we determine the Hartree-Fock approximation for the single-particle Green function in the "box- 
basis." For that purpose the Hartree-Fock Hamiltonian is diagonalized 

E (^^1 ^-^^ + ^r Ifc») (fc"l")y = ^Sf (fc.l«>y ■ (18) 

n— 1 

Here and in the following the set of basis states in the box has been truncated by assuming an appropriate -/V,nax- In 
the basis of Hartree-Fock states |a), the Hartree-Fock propagator is diagonal and given by 



-,(0)/ 



1 



where the sign in front of the infinitesimal imaginary quantity irj is positive (negative) if e^j^ is above (below) the 
Fermi energy. With these ingredients one can solve the Dyson equation (m. One possibility is to determine first the 
so-called reducible self-energy, originating from an iteration of AE, by solving 

(a| J::f{E) 1/3) = {a\ AI]y(i?) |/3) + ^ H ASy (i?) |7> ft? (7; E) (7I ^^l^'iE) |/3) (20) 
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and obtain the propagator from 

gi,{a,P;E)=gl°\a;E){a\J:^;'{E)\f3)gl°\P;E). (21) 

Using this representation of the Green function one can calculate the spectral function in the "box basis" from 

Sf^{h,E) = ;^Ini I }_J (fc,|a)y gi,{a, /3; E) {/3\h),^ ) . (22) 

For energies E below the lowest single-particle energy of a given Hartree-Fock state (with Ij) this spectral function is 
different from zero only due to the imaginary part in Y,^^'^. This contribution involves the coupling to the continuum 
of 2hlp states and is therefore non-vanishing only for energies at which the corresponding irreducible self-energy 
AE has a non-zero imaginary part. Besides this continuum contribution, the hole spectral function also receives 
contributions from the quasihole states [pi . The energies and wavefunctions of these quasihole states can be determined 
by diagonalizing the Hartree-Fock Hamiltonian plus AE in the "box basis" 

E (fc'^l ^'5,„ + J:^/ + AEy(£; = e%) \K) (fc„|T),^. = e% {h\T),^ .. (23) 

n=l 



Since in the present work AS only contains a sizeable imaginary part for energies E below e^ , the energies of the 
quasihole states come out real and the continuum contribution to the spectral function is separated in energy from 
the quasihole contribution. The quasihole contribution to the hole spectral function is given by 



with the spectroscopic factor for the quasihole state given by g 



'■^Tlj 



^ d{T\/^^i,{E)\T) 



dE 



-rij 



(24) 



(25) 



Finally, the continuum contibution of Eq. (|2^) and the quasihole parts of Eq. (^^ , which are obtained in the basis 
of box states, can be added and renormalized to obtain the spectral function in the continuum representation at the 
momenta defined by Eq. ([15h 
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Sij ih,E)^--^ {Sf^ {h ,E)+^ S% {h , E)) 



ttN, 



(26) 



C. Ground-state Properties 



The single-particle propagator calculated by the techniques described above, may also be used to evaluate expec- 
tation values of single-particle operators, like the mean square radius, and the energy of the ground state. For that 
purpose one also needs the non-diagonal part of the density matrix, which is given in the "box basis," defined in the 
preceeding subsection, by 



nij [k 



,fc„)= r dE-ImlY,{h 

"" \a.j3 



0')ij9ij{oi,P;E) {P\K)ij , 



(27) 



and contains, as before in the case of the spectral function, a continuous contribution and a part originating from the 
quasihole states 



~ qh 



{ki, fc„) = ^ Zrij {ki\T)i^ (T|fc„ 



'ij 



(28) 



The sum in this equation is restricted to quasihole states with energies below the Fermi energy ep. With this density 
matrix the expectation value for the square of the radius can be calculated according to 



Afir 



«|r2|0=^2(2j + l) J2 {h\r^\kn)inij{h,k^) 

1,3 



(29) 



i,n— 1 



with a factor of 2 accounting for isospin degeneracy. The matrix elements for r^ are given by 



(fc^l r' \k„)i = NuNru / drr^ji{hr)ji{knr) 
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(30) 



In the same way one can also calculate the expectation value for the particle number. The total energy of the ground 
state is obtained from the "Koltun sum rule" 



ii;o^ = ^2(2j + l) J2 r '^A (S + ^) (^y(fc-^) + E^Ty(fc.'^)) 



(31) 



As in Eq. (26), the sum over quasihole states T is restricted to those below ep- 



III. RESULTS AND DISCUSSION 

In our discussion of the hole spectral function in the preceeding section we have distinguished the contributions 
originating from the quasihole states and the continuum of 2hlp configurations (see Eq. (pq)). This separation 
into the two parts is displayed in Fig. y for the energy integrated spectral function (including all energies below the 
Fermi energy ep) for different partial waves Ij. This figure displays quite clearly that the momentum distribution at 
small momenta is dominated by the quasihole contribution (for those partial waves for which it exists) whereas the 



eaving the remaining nucleus (e.g. N 



high- momentum components are given by the continuum part (see also Ref. |26 ). 

This implies that a nucleon knockout reaction with small energy transfer, 
in the present case as all results presented here refer to "'^^O) in its ground state or in the lowest state with angular 
momentum and parity defined by the partial wave quantum numbers j and /, should display a spectral distribution as 
presented by the quasihole part. The high-momentum components of the spectral function (or momentum distribution) 
should only be observed in experiments which also include knockout processes into states represented by the 2hlp 
continuum. We recall that the present approach has been designed to account for the effects of short-range correlations. 
Effects due to configuration mixing of the hole state with 2hlp configuration at low energies must be treated in terms 
of shell- model configuration mixing or by techniques as discussed in Refs. [p9|-p2|. 

In order to characterize the energy dependence of the spectral functions one may define a mean value for the energy 
of the 2hlp continuum for each momentum and each partial wave by 

fdEESfJk.E) 

Typical values for this mean value range from -80 to -150 MeV for the momenta k considered in this analysis (fc < 3.3 
fm~^). One also finds that this mean value is quite independent of the partial wave considered (see left part of Fig. 
|3|) . Therefore it is useful to define a mean value of the energy by averaging over all partial waves 

^^" Eui^J + '^) IdESf^{k,E) ■ ^'^^ 

The resulting energy spectrum £{k) is shown in the left part of Fig. and compared to a simple parametrisation of 
this curve in terms of —k'^/{2m*) — C with m* = 2400 MeV and C — 80 MeV. This parametrisation demonstrates 
that the momentum dependence of this mean value is weak as compared e.g. to the kinetic energy. One may also 
compare the mean value £{k) determined by Eq. ( p3| ) in ^®0 with the corresponding quantity obtained for nuclear 
matter using the Reid potential |ljJ34]. The mean value calculated for nuclear matter shows a stronger momentum 
dependence and therefore, at high momenta, yields energies considerably below those displayed in Fig. y. This implies 
that the nuclear matter calculation exhibits a larger probability to excite 2plh configurations at higher energies as 
compared to the present approach. We will come back to a discussion of possible differences between the present 
calculation and studies in nuclear matter when we analyze the results for the momentum distribution below. 

In order to show the importance of the continuum part of the spectral functions as compared to the quasihole 
contribution and to visualize the effects of correlations, we have included in Tabfl the particle numbers for each 
partial wave including the degeneracy of the states 

dE I dkk^Sij{k,E) , (34) 

-oo Jo 

also separating the contributions originating from the quasihole states and those due to the continuum. 

In the present approach the quasihole states, which in a Hartree-Fock approximation would be occupied with a 
probability of 1.0, are occupied with a probability of 0.78, 0.91 and 0.90 in the case of Si /2, P3/2 and Pi/2, respectively. 
This means that only 14.025 out of the 16 nucleons of ^®0 occupy the quasihole states. Another 1.13 "nucleons" 
are found in the 2hlp continuum with partial wave quantum numbers of the s and p shell, while an additional 0.687 
"nucleons" are obtained from the continuum with orbital quantum numbers of the d and / shells. The distinction 
between quasihole and continuum contributions is somewhat artificial for the S2/2 orbital since the coupling to low- 
lying 2hlp states leads to a strong fragmentation of the strength pS] , which is also observed experimentally Pq| . A 
recent (e, e'p) experiment on ^^O ||37| has provided detailed information on the spectroscopic factors at low-energy 
transfer. The analysis of the experiment indicates that e.g. the Pi/2 quasihole state carries only 63% of the strength. 
This result should be compared to the 90% obtained here. This discrepancy is partly due to the emphasis in the present 
work on the accurate treatment of short-range correlations. Long-range (low-energy) correlations, not considered in 
this work, typically yield another 10% reduction of the quasihole strength 129-^2^3]. It has also been observed that a 



correct treatment of the center of mass motion may be responsible for another 10% reduction in the quasihole strength 

11- 

The sum of the particle numbers listed in Tabflis slightly smaller (15.841) than the particle number corresponding 
to ^^O. There are several possible sources for this discrepancy: First of all our analysis only accounts for momenta k 
below 3.3 fm~^ and we did not consider partial waves with I > 3. The restriction in k is determined by the choice 
of iVinax in truncating the "box basis" (see e.g. Eq. (|lq)). Inspecting the decrease of the occupation numbers listed 
in Tab.^ with increasing I one can expect that the "missing" nucleons may be found in partial waves with / > 3. 
Furthermore, however, one must keep in mind that the present approach to the single-particle Green function is not 
number-conserving, as the Green functions used to evaluate the self-energy are not determined in a self-consistent way 
0. It should be pointed out that the depletion of the occupation probabilities of the hole states, indicated in Tab.|, 
is particularly large for the Si/2 orbit. This feature can be ascribed to the closeness of the Si/2 Hartree-Fock energy 
to the 2hlp continuum which yields more leakage of strength to the continuum than for the pi/2 and P3/2 quasihole 
states. 

Inspecting the contributions to nf ■ originating from the various energy regions in Tab.fl one can sec that the major 
contributions are obtained from energies around -100 MeV. Only small contributions come from energies below -150 
MeV. The same feature is also obtained if one analyzes the momentum integrated spectral function of the continuum 

/•oo 

N-{E) = V 2(2j + 1) / dk k^ S'i^{k, E) , (35) 

shown in Fig.|[ As a function of the energy of the 2hlp states, this density of states rises very rapidly just below our 
threshold for 2hlp configurations at « —40 MeV, shows a maximum at -60 MeV, a second local maximum around 
-85 MeV, reflecting possibly some shell structure, and smoothly vanishes at lower energies. This density of states 
corresponds to a prediction of the total spectral strength to be observed in knockout reactions as a function of the 
energy transfer. 

The contribution of the 2hlp continuum to the momentum distribution is presented in Fig. |5|, exhibiting the 
contributions from partial waves with various I. The momentum distributions displayed in this figure contain the 
degeneracy factors 2(2j -f 1) and are normalized in such a way that J dkni{k) yields the total number of particles 
with orbital angular momentum I in the 2plh continuum. This figure also shows that the largest contributions are 
obtained for I — 0, although the degeneracy factor is small in this case. One can see, however, that the decrease of 
the contributions with increasing I is slow, supporting the above argument that the missing particle number exhibited 
in Tab.|| should be obtained from partial waves with / > 3. In addition, the centroid of the momentum distribution is 
shifted to higher momenta with increasing I. At momenta fc w 3 fm^^ the largest contribution is obtained from / = 3. 

The total momentum distribution, including the contribution from the quasihole states is shown in Fig. |6|. This 
distribution is presented for various energy cut-offs. The quasihole part reflects the cross section for knockout reactions 
with small energy transfer, i.e. leading to the ground state of the flnal nucleus and excited states up to « 20 MeV. 
The curve denoted hy E > -100 MeV should reflect the momentum distribution including all states of the flnal nucleus 
up to around 80 MeV, etc. As has been discussed already in connection with the spectral functions of Fig. (see 
also Ref. ||26|] ), the high- momentum components of the momentum distribution due to short-range correlations are 
expected to be observable mainly in knockout experiments with an energy transfer of the order of 100 MeV. 

The total momentum distributions resulting from the quasihole states and the 2hlp continuum are displayed again 
in Fig. and compared to predictions from studies in nuclear matter |lj,p41 . In order to enable the comparison with 
the nuclear matter results, the momentum distributions resulting from the present studies have been divided by the 
particle number and are normalized in this flgure such that J d^kn{k) yields 1. This comparison demonstrates that 
the enhancement of the momentum distribution predicted by the present study for high momenta is well below the 
corresponding prediction derived from nuclear matter. 

At flrst sight this discrepancy seems to be in contradiction to the success of the Local Density Approximation found 
in Ref. ||22l. Before we reach this conclusion, however, one must consider the following points: (i) Unfortunately, 
we cannot compare the momentum distribution obtained in our study of ^®0 using the OBE potential B of p4| ] 
with a momentum distribution derived for the same interaction in nuclear matter. The curve displayed in Fig. has 
been evaluated for the Reid soft-core potential. The modern OBE potentials are considered to be "softer" than the 
older Reid potential. Therefore part of the discrepancy might be explained by the different interaction. It would be 
very useful to pursue whether various realistic interactions, indeed, predict differences in the momentum distribution, 
which might be observed in experiment, (ii) The momentum distribution of nuclear matter has been evaluated for 
the empirical saturation density. In order to compare with a momentum distribution of a light nucleus, like ^^O, 
the momentum distribution of nuclear matter at around half the saturation density would be more appropriate. 
The momentum distribution of nuclear matter tends to be smaller at high-momentum transfers for smaller densities 
|p8[ . (iii) In our present study of finite nuclei we only consider contributions to the self-energy of the nucleons up 



to second order in the G-matrix (see Fig. ^, whereas the study in nuclear matter accounts for a self-consistent 
treatment of all ladder diagrams. It is possible that a perturbative approach underestimates the high momentum 
components in the distribution, since the G-matrix is soft as compared to the bare potential, (iv) Our present approach 
underestimates the effect of low-energy excitations (see discussion of the single-particle spectrum in calculating the 
self-energy following Eq. (w). For a finite nucleus it is quite possible that an enhancement of these correlations 
due to low-energy excitations will provide an enhancement of the momentum distribution around k = 3 fm^^. (v) 
Finally, we would like to recall that partial waves with / > 3, which were ignored in the present study may provide a 
non-negligible contribution to the momentum distribution at high momenta (see also Fig. |g). 

Finally, we would like to discuss the effects of correlations which are taken into account in thepresent investigation 
beyond the BHF approximation, on the ground-state properties of ^^O. For that purpose Tab.|l| lists the ingredients 
for calculating the total energy of the ground state according to Eq. (pl|). Furthermore we present results obtained 
for the radius of the nucleon distribution (see Eq. (p9|)). 

As a first approximation we consider the Hartree-Fock (HF) approximation, which means that the self-energy of 
the nucleons is approximated by Eq. (ph. This implies that the occupation probabilities are equal to 1 for the three 
hole states Si/2, P3/2, P1/2 and otherwise. The resulting binding energy per nucleon (-1.93 MeV) is quite small. We 
believe that this small binding energy is due to the use of the nuclear matter G-matrix calculated at the saturation 
density, which overestimates the Pauli effects as compared to a BHF calculation directly for ^^O. 

The treatment of the Pauli operator is improved by adding the 2plh part (Eq. uu) minus the correction term of 
Eq. ([ll| ) to the self-energy, an approximation which we will call Brueckner-Hartree-Fock (BHF) in the following. 
Note that the occupation probabilities of the BHF approach are identical to those of HF. Indeed, this correction 
increases the calculated binding energy to -4.01 MeV. This number is in reasonable agreement with self-consistent 
BHF calculations performed for ^^O using the same interaction [p9l. However, as the single-particle states of BHF 
are more bound than the single-particle states obtained in HF, the gain in the binding energy from HF to BHF is 
accompanied by a reduction of the calculated radius of the nucleon distribution. This is the well-known phenomenon 
of the so-called "Coester band" in finite nuclei |^, which plagues microscopic attempts to calculate ground-state 
properties of nuclear systems already for a very long time [ [40[ . 

The inclusion of the 2hlp contributions to the self-energy in the complete calculation reduces the absolute values 
of the quasihole energies (compare BHF and "Total" in Tab.y). This is to be expected from our discussion following 
Eq. (13). Despite this reduction of the quasihole energies, however, the total binding energy is increased as compared 
to BHF. This increase of the binding energy is mainly due to the continuum part of the spectral function. Comparing 
the various parts of the "Koltun sumrule" of Eq. (|3^) one finds that only 37 percent of the total energy is due to 
the quasihole part of Eq. (|3l|). The dominating part (63 percent) results from the continuum part of the spectral 
functions although this continuum part only "represents 11 percent of the nucleons" (see Tab.B). 

The calculation of the radius, however, is dominated by the quasihole contribution to the density. As the quasihole 
terms have reduced energies as compared to BHF, it is plausible that the calculated radius increases in the total 
calculation as compared to BHF. Therefore the inclusion of 2hlp terms increases the calculated binding energy and 
radius, moving the results for the ground state off the Coester band into the direction of the experimental data. 
This effect is large enough to explain the discrepancy obtained between the experimental data and the results of 
microscopic Dirac-Brueckner-Hartree-Fock calculations including relativistic effects |4^. We note that inclusion of 
three-body forces in variational calculations for '^^O also yields very good results for the binding energy p2[ . 

IV. CONCLUSIONS 

An attempt has been made to derive the spectral function and the momentum distribution from a realistic OBE 
interaction directly for a finite nucleus without the assumption of a local density approximation. The correlations 
taken into account beyond the Hartree-Fock approximation yield a strong enhancement of the momentum distribution 
at high momenta. It is demonstrated that this enhancement originates from the spectral function at large negative 
energies and therefore should be observed in nucleon knockout reactions with large energy transfer leaving the final 
nucleus at an excitation energy of about 100 MeV. 

The enhancement of the high- momentum components is weaker as obtained in studies of nuclear matter. This 
difference may be due to the different interactions employed (unfortunately no nuclear matter result is available for 
the OBE interaction used here) or due to approximations used in the calculation for the finite system. Therefore 
further studies of these approximations (poor treatment of low-energy excitations, the self-energy of the nucleons is 
calculated in a perturbative way including terms up to second order in G) is required before conclusions about the 
validity of the Local Density Approximation relating the results of nuclear matter to those of finite nuclei can be 
drawn. 



The resulting Green function is also used to determine the total energy and the radius of the nucleon distribution. 
It is demonstrated that the inclusion of 2-hole 1-particle contributions to the self-energy of the nucleon yields an 
enhancement of the calculated binding energy per nucleon (w 1 MeV) and an increase of the radius (w 0.05 fm) 
for ^^O as compared to the Brueckner-Hartree-Fock approach. This could be sufficient to explain the discrepancy 
remaining between experimental data and microscopic Dirac-BHF calculations pl[ . 
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TABLE I. Distribution of nucleons in ^^O. Listed are the total occupation number h for various partial waves (see Eq.(34)) 
but also the contributions from the quasihole (n''') and the continuum part iff) of the spectral function, separately. The 
continuum part is split further into contributions originating from energies E below -150 MeV {n'^{E < —150)) and from 
energies below -100 MeV. The last line shows the sum of particle numbers for all partial waves listed. 
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TABLE II. Groundstate properties of ^®0. Listed are the energies e and kinetic energies t of the quasihole states (qh) and 
the corresponding mean values for the continuum contribution (c), normalized to 1, for the various partial waves. Multiplying 
the sum: l/2(t -I- e) of these mean values with the corresponding particle numbers of Tab. I, one obtains the contribution AE 
to the energy of the ground state (see the Koltun sumrule Eq. (31)). Summing up all these contributions and dividing by the 
nucleon number yields the energy per nucloon E/A. Furthermore we give the radius for nucleon distribution (r), calculated 
from the square root of Eq. (29). Results are presented for the Hartree-Fock (HF), Brueckner-Hartree-Fock (BHF) and the 
complete calculation (Total). The particle numbers for the qh states in HF and BHF are equal to the degeneracy of the states, 
all other occupation numbers are zero. The results for the radii are given in fm, all other entries in MeV 
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FIG. 1. Graphical representation of the Hartree-Fock (a), the 2-particle 1-holc (2plh, b) and the 2- hole 1-particle contribution 
(2hlp, c) to the self-energy of the nucleon 

FIG. 2. Momentum distribution for different partial waves in ^®0 (see Eq. (2)) obtained by integrating the spectral function 
Sijik, E) (see Eq. (3)). The momentum distribution is the sum of the quasihole contribution (dashed curve) and the continuum 
contribution (dotted curve). All functions are normalized such that ( dk n{k) = 1 if S{k) refers to an orbit which is mostly 
occupied. 

FIG. 3. Mean value for the energy of the 2hlp continuum as a function of the momentum k. The left part of this figure shows 
results for £ij (see Eq. (32)) in various partial waves. In the right part of the figure the mean value averaged for the various I 
and j (see Eq. (33)) is displayed. For comparison this part also includes a simple parametrization in terms of — fc^/(2m*) — C 
with m* = 2400 MeV and C = 80 MeV. 

FIG. 4. Density of states or total occupation probability of the 2hlp continuum as a function of the energy E (see Eq. (35)). 
The normalization of this distribution is such that the integration over the energy yields the total particle number of 1.816 (see 
Tab. I) in the continuum. 

FIG. 5. The momentum distribution for various orbital angular momenta. These distribution account for the different j, 
include the degeneracy factors 2(2j + 1), and are normalized in such a way that J dkniik) yields the total number of particles 
with orbital angular momentum I in the 2plh continuum. 

FIG. 6. The total momentum distribution of ^®0. Shown are also the quasihole contribution and the results obtained with 
various energy cut-offs in the integration of the spectral functions 

FIG. 7. The total momentum distribution obtained in the present investigation for ^®0 employing the OBE potential B of 
Ref. 24 is compared to the momentum distribution obtained for the Reid soft-core potential in nuclear matter [14]. In this 
figure the momentum distributions are normalized in such a way that f cPkn{k) yields 1. 
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